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In the present paper, we investigate the efficacy of Dzyaloshisnhkii-Moriya (DM) interaction to 
convert the bound entangled states into free entangled states. We consider the tripartite hybrid 
system as a pair of non interacting two qutrits initially prepared in bound entangled states and one 
auxiliary qubit. Here we consider two types of bound entangled states investigated by Horodecki. 

The auxiliary qubit interacts with any one of the qutrit of the pair through DM interaction. We 
show that by tuning the probability amplitude of auxiliary qubit and DM interaction strength one 
can free the bound entangled states, which can be further distilled. We use the reduction criterion 
to find the range of the parameters of probability amplitude of auxiliary qubit and DM interaction 
strength, for which the states are distillable. The realignment criterion and negativity have been 
used for detection and quantification of entanglement. 

I. INTRODUCTION 

Quantum entanglement mm is a physical phenomenon which takes place when particles interact in microscopic 
world in such a way that the quantum state of one particle can be described in terms of each other. The entangle¬ 
ment phenomenon is expected to be a useful resource for future quantum technologies. Many applications, based on 
entanglement, have been investigated, like quantum teleportation [3], quantum imaging [4], quantum game theory [5], 
secure key quantum transmission [^, etc. To develop quantum technologies, based on entanglement, the quantum 
community needs long time entangled quantum systems, which are free from noise. Quantum systems are too evasive 
as they may loose their entanglement by external environmental interactions and can go for entanglement sudden 
death (ESD) for finite time [ll[8]. So, dynamics of entanglement and its control under various environmental interac¬ 
tions conceptually underpinning the quantum information processing. Dynamics in various quantum spin chains have 
been studied under Dzyaloshisnhkii-Moriya (DM) interaction [TTJ[T2]. DM interaction is a useful resource in quantum 
information processing to entangle and disentangle the quantum systems. Recently, Zang et al. studied the entangle¬ 
ment dynamics of two qubit pair by taking an third qubit or qutrit, which interact with a qubit of the pair through 
DM interaction [IM3- They studied the dynamics of entanglement by taking a third qubit as a controller qubit. 
Eurther, they studied the same by taking a third qutrit as a controller qutrit. They proposed that by manipulating the 
probability amplitude of third qubit or qutrit and DM interaction strength one can control the entanglement between 
two qubits induced by DM interaction and hence in various quantum spin chains |16j . At this point, we mention here 
that the same method may not only be used to control the entanglement, but it may also be used to free the bound 
entanglement [nMO] in various bound entangled states. Once the bound entangled states are free, then they can be 
distilled [21]. 

The quantum states beyond the dimension 3 0 3 have been classified in two categories like free entangled states 
and bound entangled states. Eree entangled states are distillable states or in other words noise free states. On the 
other hand, bound entangled states are noisy states and no pure entanglement can be obtained by local operations 
and classical communication. It is difficult to use bound entangled states directly for practical quantum information 
processing. However, by providing additional resource, bound entanglement can be activated to increase the fidelity 
of quantum teleportation pMMj . Up to now we don’t have satisfactory tools to quantify and detecting the bound 
entanglement. Recently, the free entanglement production from bound entangled states is proposed by using the 
ancillary system which is coupled to the initial bound-entangled state via appropriate weak measurements [25] . 

In the present work, we consider qutrit-qutrit bound entangled bipartite states proposed by Horodecki UMOj- The 
main goal of the present work is to show that DM interaction can be a used as a useful agent to free the bound 
entangled state in two qutrit system. Once the states are free, further these can be distilled. We have used reduction 
criterion to check the distillability of these states. Using this criterion, one can ensure that the states are surely 
distillable. Further, we have used realignment criterion to detect the bound entanglement and negativity to measure 
the free entanglement. The motivation of this study comes from our recent works P51I29| . 

The plan of the paper is as follows. In Sect. 2, we present the Hamiltonian of the system. Sect 3. is devoted to the 
description of Horodeckis bound entangled states and reduced system dynamics. In Sect. 4, we discuss the reduction 
criterion, realignment criterion and negativity. Sect 5. is devoted to time evolution of negativity and realignment 
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FIG. 1: Plot of entanglement and realignment criterion with D = 0.0. Green color graph represents negativity (N) and red 
color graph represents the realignement criterion (R). 


criterion. Lastly in Sect. 6, we report our conclusion. 


II. HAMILTONIAN OF THE SYSTEM 

We consider a qutrit (A)-qutrit (B) pair and an auxiliary qubit (C) which interact with any one of the qutrit of the 
pair through DM interaction. Here, we assume that the auxiliary qubit (C) interact with the qutrit (B) of the pair. 
Now the Hamiltonian of the system can be written as 


( 1 ) 

where Hab is the Hamiltonian of qutrit (A) and qutrit (B) and H'‘gQ is the interaction Hamiltonian of qutrit (B) and 
qubit (C). We consider uncoupled qutrit (A) and qutrit (B), so Hab is zero. Now the Hamiltonian becomes 

H = H^^l. = D.{aBXac), ( 2 ) 

where D is DM interaction between qutrit (B) and qubit (C) and ah, dc are associated vectors of qutrit (B) and 
qubit (C) respectively. We assume that DM interaction exist along the z-direction only. In this case, the Hamiltonian 
can be simplihed as 


H = D.{ag ® a^ — a]^ ® a^), (3) 

where tr^ and a^ are Gell-Mann matrices for qutrit (B) and a^ and cr^ are X and Y Pauli matrices of qubit (C) 
respectively. The above Hamiltonian is a matrix having 6x6 dimension. Further, it is multiplied by the identity 
matrix of dimension three and we obtained the dimension as 18 x 18. The matrix of Hamiltonian is easy to diagonalize 
by using the method of eigendecomposition. The unitary time evolution operator is easily commutable as 

U{t)=e-^^\ (4) 


which is also a 18 x 18 matrix. This matrix has been used to obtain the time evolution of density matrix of the system. 


HI. HORODECKI’S BOUND ENTANGLED STATES AND REDUCED SYSTEM DYNAMICS 

In this section, we describe two known bipartite Horodecki’s bound entangled states niHnii] in 3 G 3 dimension. In 
subsection 3.1, we present the state 1 and in subsection 3.2 the state 2 is presented. 
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FIG. 2: Plot of negativity (N) and realignment criterion (R) with the parameter values cq = 0 and Dt = 0. Green color graph 
represents negativity (N) and red graph represents realignment criterion (R). 


A. State 1. 

The Horodecki’s bound entangled state [HI HI] in 3 0 3 dimension is given by 

Pa{Q) = jP + + 21 2 < a < 5, 

where 

P=mw, |^) = ^(|00) + |ll) + |22)) 

Q=(|01)(01| + |12)(12| + |20)(20|), 
i?= (|10)(10| + |21)(21| + |02)(02|). 

Horodecki demonstrated that 


( 5 ) 

( 6 ) 

( 7 ) 

( 8 ) 


{ Separable for 2 < a < 3, 

Bound entangled for 3 < a < 4, (9) 

Free entangled for 4 < a < 5. 


B. State 2. 


Another well known bound entangled state investigated by Horodecki [MUD] in 3 0 3 dimension is given as 


Pa 


1 

8a + 1 


a 0 0 0 a 0 
0 a 0 0 0 0 
0 0 a 0 0 0 
0 0 0 a 0 0 
a 0 0 0 a 0 
0 0 0 0 0 a 
0 0 0 0 0 0 
0 0 0 0 0 0 
a 0 0 0 a 0 


0 0a 

0 0 0 

0 0 0 

0 0 0 

0 0a 

0 0 0 

l+a n Vl — 

2^2 
0 a 0 

Vl— n 1+g 
2^2 


( 10 ) 
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here 0 < a < 1. 


C. Reduced system dynamics 

In this subsection, we obtain reduced density matrix of qutrit-qutrit system. To begin with, let us consider the 
auxiliary qubit (C) in pure state as 


k) = co|0)+ci|l) (11) 

with normalization condition 

|coP + |ciP = l (12) 

where the probability amplitudes cq and ci are complex in general. 

Qutrit (A)-qutrit (B) is prepared initially in bound entangled state. The density matrix of bound entangled state 
before interaction with qubit (C) is given by ([^. Now we can obtain the composite density matrix of the open system, 
Pcomp.i^) 


Pcomp.(O) — /^a(0) ^ Pc: (^^) 

where pc is the density matrix of the auxiliary qubit (C). The density matrix after the interaction at time t is given 

by 


PcompX^) — ^XXP comp. mHt), (14) 

where U(t) is unitary time evolution operator given by Now we make the order same of all the matrices involved 
in equation (14) as 18 x 18. This can be done by doing the tensor product of Pcomp.(O) with identity matrix of the 
order 2x2. Now the order of the matrix Pcom.{t) becomes as 18 x 18. Taking partial trace of Pcom.X) over the basis 
of auxiliary qubit (C), we get the reduced density matrix of 9 x 9 dimension as, 

= Ptrc[pco7np.{t)]. (15) 

Now we obtain the reduced density matrix for bound entangled state proposed by the Horodecki given by Eq. ([^. 
The reduced density matrix is given below 
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FIG. 3: Plot of eigenvalues (Eig.) of reduction criterion with the parameters range 0 < Dt < 5 and 0 < a < 1. 


where 
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(16) 


and 

p = sin[-\/2ZIf]coCi, q = cos[-\/2£)f], r = ^ sin[2-\/2IIt]coCi, s = ^ cos[2y/2Dt]coCi. 

Observing the reduced density matrix, we conclude that cq, D and t are the key parameters which influence the 
entanglement between two qutrit pair. Similarly, we obtain the reduced density matrix corresponding to another 
bound entangled state given in subsection IIIB with the help of Eq. (15). 


IV. REDUCTION, REALIGNMENT CRITERION AND NEGATIVITY 

In this section, we discuss the reduction criterion and its ability to detect the separability and distillability of 
bipartite states. For the separability of a composite system AB, the following conditions should be satisfied 


PA®I — p>^ and I^pB—p>0, 


( 17 ) 
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FIG. 4: Contour plot of Fig 2, the state is distillable with the parameters range 1.59 < Dt < and 0 < q < 3.439 


where pA =TrB{p), Pb = TrA{p) and p is the composite density matrix of the system AB. The states which violates 
the condition 0 can be distilled. So, for distillability the following condition should be necessarily satisfied [30| . 

PA^I — p<0 or I®pB — p<Q- (18) 

Further, we obtain the partial transpose and realigned matrix of the reduced density matrices to obtain the negativity 
and realignment by using the following matrices 


(19) 

( 20 ) 

These matrices have been used to calculate the quantities defined as 


{Pij,kl) ~ Pil,kj 

p R \ 

\Pij,kl) ~ Pikjl- 


N = 


IIP^I|-1 

2 



( 21 ) 


where 11.. 11 is the trace norm of the matrix. The first quantity corresponds to negativity |31H35) and it has been 
used to measure the free entanglement while the second one has been used to detect the bound entanglement in the 
qutrit-qutrit system. Either > 0 or i? > 0 implies that the state is entangled. N = 0 and R > 0 implies that the 
state is PPT bound entangled, and N > 0 corresponds to free entangled state. However, there is no evidence for the 
existence of NPT bound entangled states [3S] as yet, but we can not avoid this future possibility. 


V. TIME EVOLUTION OF NEGATIVITY AND REALIGNMENT CRITERION 

In this section, first we verify the calculation of reduced density matrices for {Dt = 0) for state 1. The reduced 
density matrix of two qutrits corresponding to state 1 is given by Eq. ( |I^ . We put the values of the parameters Dt 
and Co as zero in reduced density matrix then reduced density matrix maps to the initial state of two qutrits given by 
Eq. (§, which verify the correctness of reduced density matrix obtained in calculation. Further, we plot the negativity 
and realignment criterion and the result is shown in Fig. ^ The green color graph represent the negativity (N), while 
the red color graph represents the realignment criteria (R) in the figure. The realignment criterion for 2 < a < 3 is 
negative, so it shows that the state is separable. For 3 < a < 4 the realignment criterion is positive, which show that 
the state is bound entangled. Further, for 4 < a < 5 the negativity graph is positive, which shows that the state is 
free entangled. These results match with the result given by Eq. (|^. 

Similarly, we verify the reduced density matrix obtained for the state 2 by putting the parameter values Dt = 0 and 
Co = 0. The reduced density matrix maps to the initial density matrix of the state 2 given in subsection [III B[ which 
proves that the reduced density matrix obtained is correct. Now, we plot the negativity and realignment criterion 
both in the absence of DM interaction in Fig. The realignment criterion achieves positive values for 0 < a < 1 
and negativity is zero. These results show that initially the state 2 of two qutrits is bound entangled with 0 < a < 1. 
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FIG. 5: Plot of negativity (N) and realignment criterion (R) with the parameter range co = 0.2 and (1.5 < Dt < 3.7, 
0 < a < 0.5). Green color graphs represent N and red graphs represent R. Realignment criterion fails at a = 0.02. 


Further, in this section we detect the distillability of the Horodecki bound entangled states 1 and 2, by using the 
reduction criterion with the reduced density matrices of two qntrits in the presence of DM interaction. Next, we study 
the evolution of free entanglement and realignment criterion obtained by using Eq. (21) for the bound entangled states 
1 and 2 in these two cases. In case 1, we consider the bound entangled state described under the subsection 3.1 and 
in case 2, we do the same for another bound entangled state given in subsection 3.2. 


A. Case 1 


In this case, we consider the state 1 given in subsection 3.1. Here, we replace c{ in Eq. (21) in terms of Cg by using 
the normalization condition given in Eq. ( |I^ . Next we check the distillability of the state by using the reduction 
criterion give in Eq. ( [l^ . We found that the matrices in reduction criterion, i.e., {pA ® I — p) or {I ® Pb — p), 
incorporate the parameters Dt, cg and a. Further, the eigenvalues of these matrices have been calculated by taking 
the range of the parameters as (3 < a < 4, 0 < Dt < 5) with the varying value of the parameter cg. With the range 
(0 < Cg < 1). It is found that the eigenvalues of these matrices achieve positive values. These eigenvalues are obtained 
with the maximum values of Cg (i.e. Cg = 1), it is found that all eigenvalues are positive, so the reduction criterion 
given in Eq. (181 fails and state 1 is not distillable. 


B. Case 2 


In this case, we consider the state 2 given in subsection 3.2. Here, we again replace cf in Eq. (21) in terms of 
Cg by using the normalization condition given in Eq. (12). Further, we check the distillability of the state by using 
the reduced density matrix corresponding to this state in the presence of DM interaction. The distillability of the 
state is detected by using reduction criterion given by Eq. (18). The matrices involved in left hand side of reduction 
criterion i.e., {pA ® I — p) or {I® pB — p), incorporate the parameters Dt, cg and a. The eigenvalues of these matrices 
have been calculated by taking the range of the parameters as (3 < a < 4, 0 < Dt < 5) with the varying values of 
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FIG. 6: Plot of negativity (N) and realignment criterion (R) with the parameter range co = 0.4 and 1.5 < Dt < 3.7. Green 
color graphs represents N and red color graph represents R. 


parameter cg in the range (0 < Cg < 1). It is found that the eigenvalues of these matrices achieve negative values and 
hence the reduction criterion given by Eq. (181 is satisfied for state 2. Our goal is to find out the maximum range 
of the parameters Dt and a for which the eigenvalues are negative. By simulation, we found these maximum ranges 
corresponding to the parameter values cg = 0.7. These eigenvalues are plotted in Fig. [^with the parameters range 
0 < Dt < 5 and 0 < a < 1.0 corresponding to the parameter value cg = 0.7. The corresponding contour plot of the 
Fig.0 is depicted in Fig. With the help of this contour plot, we obtained the maximum ranges of the parameters 
Dt and a for which the state 2 is distillable. These ranges are 1.59 < Dt < 3.75 and 0 < a < 3.439. The distillability 
of the state proves that the states are free entangled. Only free entangled state can be distilled in Bell like pairs, which 
may be used, further, for quantum information applications. So, based on the parameters ranges 1.59 < Dt < 3.755 
and 0 < Of < 3.439, we plot the negativity (N) and realignment criterion (R) obtained from Eq. (21) with different 
values of the parameter cg in Figs. and[^ First, we fix the value of the parameter cg = 0.2 and vary the value of 
Dt within the range 1.5 < Dt < 3.7 and for a with 0 < a < 3.439 and plot the results in Fig. Observing Fig. 
it is concluded that with a = 0 the realignment criterion achieves the positive values, so initially the state is bound 
entangled, but as the value of a advances the states become free. So these states can be easily distillable. As the 
value of the parameter Dt increases, the initial amplitudes of both realignment and negativity fluctuate between 0.04 
and 0.05. At a particular value of a = 0.02, the realignment criterion becomes negative, so it fails to detect the bound 
entanglement. But corresponding to negative realignment criterion, we can not avoid the possibility of NPT bound 
entangled states. We observe that for the parameter values, cg = 0.2, Dt = 3.0 the free entanglement in the states 
vanish after a = 0.12. 

Next, we increase the value of the parameter cg as 0.4 and sketch the graphs between Dt and a in Fig. We 
found that as the value of parameter cg increases from 0.2 to 0.4, the initial amplitudes of both realignment criterion 
and negativity increases and fluctuates within the limits 0.08 to 0.10. Initially the states are bound entangled but 
as the value of the parameter Dt increases the states become free and hence can be distilled. Again, we found 
that for some values of a, the realignment criterion fails to detect the bound entangled states. These values are 
{a = 0.07,0.04,0.08,0.06). After these values the realignment criterion becomes negative and the possibility of NPT 
bound entangled states can not be avoided. 

Further, we continue our study for cg = 0.7 and the results are shown in Fig. Initially the states are bound 
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FIG. 7: Plot of negativity (N) and realignment criterion (R) with the parameter range co = 0.7 and 1.5 < Dt < 3.7. Green 
color graphs represent N and red color graphs represent R. 


entangled as the realignment criterion is positive for a = 0. As the value of parameter a advances, the states become 
free and hence distillable. The amplitudes of both negativity and realignment criterion fluctuate between the range 
0.12 to 0.20. We have found the situation when the states becomes totally free and fully distillable. These parameter 
values are (cq = 0.7, Dt = 1.5), (cg = 0.7, Dt = 2.5) and (cg = 0.7, Dt = 2.8). However, we detect the values of 
the parameter a where the realignment criterion becomes negative, these values are a = 0.07,0.24,0.14. After these 
values of a, the realignment criterion fails and corresponding to negative portion of realignment criterion possibility 
of NPT bound entangled states can not be avoided. 
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VI. CONCLUSION 

In this article, we presented a method to distill the bound entangled states by using DM interaction. We consider 
two qutrits initially prepared in Horodecki bound entangled states and one auxiliary qubit. The auxiliary qubit is 
prepared in pure state, which interact with any one of the qutrit through DM interaction. By varying the probability 
amplitude of the auxiliary qubit and DM interaction strength the states can be free and further can be distilled. 
We have used the reduction criterion, which is necessary condition for distillability, later the maximum values of the 
parameters of Cq and Dt has been obtained for which the the states are distillable. The realignment criterion is used 
to detect the bound entangled nature of the states and negativity is used to measure the free entanglement in the 
states. By varying the values of the parameters cq, a and Dt, the Horodecki bound entangled states can be converted 
in free entangled states. Free entangled states are easily distillable. We have found that DM interaction can be used 
to free the Horodecki bound entangled state 2. With certain parameter values of Cq and Dt, it is completely free and 
hence distillable. These values are (cq = 0.7, Dt = 1.5), (cq = 0.7, Dt = 2.5) and (cq = 0.7, Dt = 2.8). As the value of 
the parameter cq increases the initial amplitudes of both negativity and realignment criterion increases. We also found 
that realignment criterion fails at particular values of a along with the values of cq and Dt. So, corresponding to the 
negative portion of the realignment criterion the possibility of NPT bound entangled states [3S] exists. We hope that 
this method of converting bound to free entanglement can be useful in quantum information processing and varieties 
of bound entangled states can be checked under DM interaction. 
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